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HIGHLIGHTS 


►  A  micro— macro  homogeneous  model  is  developed  for  a  rechargeable  Li— air  battery. 

►  The  model  includes  the  mass  transport  through  the  lithium  oxide  layer. 

►  The  model  can  simulate  the  discharge  and  charge  feature  without  parasitic  reactions. 

►  Effects  of  cathode  structure  and  feeding  conditions  have  been  investigated. 
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A  micro— macro  homogeneous  mathematical  model  is  developed  for  a  rechargeable  Li— air  battery  using 
a  concentrated  binary  electrolyte  theory.  The  dynamic  behaviour  of  the  porous  cathode  is  determined  by 
a  numerical  solution  of  the  combined  continuity,  transport  and  kinetics  equations.  This  model  considers 
the  time  and  space  dependence  of  the  battery  system,  including  the  mass  transport  along  the  depth  of 
the  cell,  the  local  mass  transfer  between  lithium  peroxide  (Li202)  layer  inside  the  cathode  and  active 
surface  morphology  changing  with  the  Li202  growth.  The  model  predicts  that  the  battery  capacity  and 
discharge  potential  are  sensitive  not  only  to  the  solubility  of  oxygen,  which  is  predominantly  limited  by 
depletion  of  the  oxygen  concentration,  as  well  as  to  the  cathode  porosity,  the  cathode  structure  and 
kinetic  parameters.  In  addition,  the  charging  behaviour  is  simulated  by  modelling  the  dissolution  of  solid 
Li202  product.  The  model  suggests  that  the  charging  voltage  can  be  lowered  depending  on  capability  of 
electrolyte  to  dissolve  the  Li202  discharge  products.  As  a  result,  a  simulation  of  electrode  behaviour  is 
obtained. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  development  of  high  performance  and  light  weight  energy 
storage  devices  has  recently  focused  on  the  rechargeable  lithium- 
air  batteries.  Mainly  due  to  their  high  energy  density,  theoretically 
up  to  11,400  Wh  kg-1  (which  is  about  10  times  greater  than  well- 
known  lithium-ion  battery),  such  batteries  are  now  considered  as 
one  of  the  promising  alternatives  to  lithium-ion  batteries  with 
potentially  wide  applications,  from  small  portable  electronics  to 
electric  vehicles.  The  first  Li— air  system  with  a  non-aqueous  elec¬ 
trolyte  was  presented  by  Abraham  and  Jiang  in  1996  [1],  After  their 
promising  rechargeable  ability  was  demonstrated  by  Bruce  group 
[2],  Li— air  batteries  have  attracted  much  more  attention  [2—7], 
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As  shown  in  Fig.  1,  a  Li— air  battery  contains  a  metal  lithium 
anode,  a  solid  polymer  separator  and  a  porous  carbon  or  catalyst- 
loaded  carbon  air  electrode  filled  with  an  organic  electrolyte 
comprising  a  dissolve  lithium  salt  in  an  aprotic  solvent.  The  porous 
carbon  electrode  provides  a  site  for  electrochemical  reduction  of 
oxygen  and,  during  battery  operation,  oxygen  (comes  from  the 
external  air)  is  dissolved  in  the  electrolyte,  penetrates  through  the 
pore  of  the  cathode  and  reacts  with  the  lithium  ion  (Li+)  at  the 
active  site.  It  is  important  to  understand  the  oxygen  reduction 
mechanisms,  including  intermediate  steps  (e.g.  discharge  products 
formation).  Up  to  date,  there  are  various  mechanisms  for  oxygen 
reduction  reaction  (ORR),  related  to  type  of  electrolyte,  catalyst  and 
battery  operating  conditions.  Peng  et  al.  [8]  investigated  in  situ 
spectroscopic  data  of  oxygen  reaction  in  a  non-aqueous  electrolyte 
and  found  strong  evidence  that  lithium  superoxide  (Li02)  is  indeed 
an  intermediate  species  during  oxygen  reduction  before  dispro¬ 
portion  to  final  Li202  product.  For  battery  charging,  the  reaction  is 
simply  the  oxidation  of  Li202  directly  into  oxygen  and  Li+  without 
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Fig.  1.  Schematic  diagram  computation  domain  of  a  Li-air  cell  during  discharge.  The 
inset  demonstrates  the  porous  carbon  cathode. 


passing  intermediate  Li02  route.  Based  on  literature  [1,9,10],  our 
model  assumes  the  following  battery  reactions: 

Anode 

Li ->  Li+  +  e  (E°  =  0.00  V)  (I) 

Cathode 


2 Li+  +  02  +  2e~  —  Li202  (£°  =  2.96  V)  (II) 

4Li+  +  02  +4e-  -*■  2Li20  (E°  =  2.91  V)  (III) 

Li202  is  the  main  reaction  product  considered  in  this  model 
because  it  is  insoluble  in  the  organic  electrolyte  and  could  cover  on 
the  active  surface  area  inside  the  porous  cathode  as  well  as  block 
the  pathway  for  reactive  species  (Li+  and  oxygen),  preventing 
further  reactions  inside  the  cathode  and  contributing  to  the  end  of 
cell  discharge  [9,10].  Conversion  of  Li202  to  Li20  is  not  considered 
because,  if  any,  it  only  occurs  at  a  significant  more  negative  (lower) 
potential  below  2  V,  which  is  beyond  the  potential  range  of  the 
model.  Thus,  a  micro-macro  homogeneous  mathematical  model, 
dealing  with  the  entire  electrode— electrolyte  cell  system  as  two 
continua  (one  from  the  electrode  matrix  and  the  other  from  the 
solution  filling  in  all  the  space  of  the  electrode),  is  developed  for  the 
cycling  operation  of  a  porous  cathode  of  a  rechargeable  Li— air 
battery  using  a  concentrated  binary  electrolyte  theory  [11,12],  The 
model  predicts  the  time  dependence  of  electrolyte  concentration, 
non-uniform  porosity  and  reaction  rate.  The  more  accurate  kinetic 
reaction  in  a  porous  cathode  is  established  by  including  the 
dependency  of  the  diffusion  transport  on  reactant  concentrations 
and  the  solubility  limit  of  the  product.  The  model  also  includes  the 
mass  transport  along  the  depth  of  the  cell,  the  local  mass  transfer 
between  Li202  layers  and  active  surface  morphology  changing  with 
the  Li202  formation.  Important  parameters,  such  as  porosity  and 
thickness  of  electrode,  electrolyte  transport  properties  and  surface 
activity,  are  investigated  through  this  modelling  study.  The  model 
predicts  the  direction  of  improving  the  Li— air  battery  performance 
in  term  of  the  discharge  product  formation  and  dynamic  model 
parameters. 

Different  from  previous  approaches  [13,14],  assuming  that  the 
cathode  only  contains  a  large  number  of  cylindrical  pores  without 
considering  the  macroscopic  aspects,  our  model  considers  elec¬ 
trochemical  kinetic  dependency  on  both  Li+  and  oxygen  species 


and  includes  local  mass  transport  through  the  solid  Li202  layer.  In 
addition,  variable  physical  and  chemical  properties  are  carefully 
applied  in  the  model  based  on  available  published  experiment 
data.  Therefore,  this  model  can  be  used  to  describe  the  behaviour  of 
Li— air  batteries  as  well  as  to  optimise  the  performance  of  these 
batteries. 


2.  Model  development 

2.1.  Basic  equations 


A  mathematical  description  for  a  prismatic  of  Li— air  battery  has 
been  developed  involving  conservation  of  mass  and  current, 
species  transport,  and  reaction  kinetic  in  the  cathode  and  separator 
to  clarify  the  mechanism  inside  the  cell.  The  model  used  in  this 
work  is  for  a  Li— air  cell  consisting  of  a  thin  lithium  sheet  negative 
electrode,  an  anode  protective  layer  (APL),  a  separator,  and  a  porous 
carbon  oxygen/air  cathode  filled  with  an  organic  electrolyte  as 
shown  in  Fig.  1.  Current  collectors  are  placed  at  the  back  of  each 
electrode.  A  concentrated  binary  electrolyte  theory  is  used  to 
describe  the  motion  for  of  each  species  in  the  electrolytic  solution. 
The  Li202  formation  inside  the  porous  cathode  is  presented  using 
macro-homogeneous  porous  model,  defining  the  electrode  by  its 
porosity  which  is  initially  uniform  but  changes  during  discharge. 
The  macroscopic  theory  of  porous  electrode  treats  the  solution  and 
solid  matrix  phases  as  superimposed  continuum.  Based  on  this 
approach,  a  material  balance  equation  for  specie  i  in  the  electrolyte 
solution  phase  can  be  expressed  as 


9(eCj) 

8t 


— V-JVj  +tj 


(1) 


where  c,  is  the  bulk  concentration  of  species  i  in  the  solution  phase 
which  is  averaged  over  the  volume  of  the  solution  in  the  pores,  e  is 
the  porosity  of  the  electrode  which  is  the  electrolyte  space  in  the 
solution  phase,  JV;  is  the  molar  flux  of  species  i  in  the  porous  solu¬ 
tion  averaged  over  the  cross  sectional  area  of  the  electrode,  and  r,  is 
the  volumetric  production  rate  of  species  i  from  the  solid  phase 
(electrode  material)  to  solution  phase  (electrolyte  in  the  porous) 
within  the  porous  electrode.  The  concentration  of  lithium  electro¬ 
lyte  is  the  same  as  concentration  of  Li+  due  to  the  binary  electrolyte 
assumption.  Without  convection  term,  the  diffusion  and  migration 
fluxes  equations  for  mass  transfers  for  both  Li+  and  oxygen  in  the 
cell  can  be  expressed  as 

Nu  =  -DLi,effVcLi  +  1-^  (2) 


N02  =  ~D02.eff^c02 


(3) 


where  F>u,eff  and  D0ieff  are  the  effective  diffusion  coefficient  of  Li+ 
and  02,  respectively,  t+  is  the  transference  number  of  Li+,  F  is 
Faraday’s  constant  which  is  equal  to  96,485  C  mol-1,  and  i2  is  the 
current  density  in  the  solution  phase  or  electrolyte  current  density 
which  can  be  defined  by  the  gradient  of  the  potential  in  concen¬ 
trated  electrolyte  solution  as  [15,16] 


V</-2  =  -  '2 
Keff 


vRT  (  s+ 
f  \nv+ 


4 

z+v+ 


nc0J\  dine ) 


(4) 


where  n,  s+,  So,  v+,  and  v  represent  the  number  of  electrons  trans¬ 
ferred,  the  stoichiometric  coefficients  for  cation  and  solvent,  the 
numbers  of  cation,  and  the  number  of  moles  of  ions  into  which 
a  mole  of  electrolyte  dissociates,  respectively,  c  and  Co  are  the  molar 
concentrations  of  the  electrolyte  and  solvent  in  the  electrolyte 
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phase,  respectively.  In  Li— air  battery,  a  1:1  binary  electrolyte  is 
applied  in  the  cell.  Then,  so  =  0,  s+  =  -1,  v+  =  n  =  1,  and  v  =  2.  The 
Eq.  (4)  can  be  rearranged  into  the  current  density  in  electrolyte  as 


h  =  -KeffVfa  -  - 


4.*q 

'dlncLiJ 


(5) 


where  Kejj  is  the  effective  conductivity  of  the  electrolyte,  <p2  is  the 
electrolyte  potential  (electric  potential  of  Li+),  R  is  the  universal  gas 
constant  which  is  equal  to  8.3143  J  mol-1  K-1,  T  is  the  cell 
temperature  in  Kelvin,  and  /is  the  activity  coefficient  of  LiPF6  salt. 

In  the  solid  matrix  phase,  the  movement  of  electron  is  governed 
by  Ohm’s  law  which  computes  the  electric  potential  or  potential  of 
electron,  <pt 

H  =  -  VeJfVfo  (6) 

where  rreff  is  the  effective  conductivity  of  the  electron  in  the  elec¬ 
trode.  This  parameter  is  affected  by  the  volume  fraction  of  solid 
electrode.  The  effective  parameters  Dejf,  Keff  and  aejj  in  the  above 
equations  are  usually  depend  on  the  tortuosity  of  individual  phase 
in  porous  cathode  (through  porosity  or  volume  fraction).  These 
parameters  are  applied  for  only  the  porous  cathode  region  (others 
are  not  porosity)  and  are  corrected  to  account  for  the  porosity  effect 
using  Bruggeman  correlation  [17] 


density.  The  value  of  s;,  z*.  and  n  can  be  defined  by  matching  with  an 
individual  electrode  reaction  using  general  form  of  Eq.  (12),  for 
example,  the  value  of  Su,  Zu,  and  n  of  Li+  from  Eq.  II  are  -2, 1,  and  2, 
respectively. 

Only  one  electrochemical  reaction,  i.e.  Li202  formation  inside 
porous  electrode,  is  considered  in  the  present  work.  In  practical 
cell,  the  battery  side  reactions  during  cycling  can  occur  along  with 
the  main  electrode  reactions.  This  is  reflected  by  the  superficial 
production  rate  of  a  species  from  solid  phase  to  pore  solution,  given 
by  Faraday’s  law 


(14) 


2.3.  Rate  expressions 

For  electrochemical  reaction  at  the  cathode,  the  Butler-Volmer 
equation  is  applied  in  the  model  using  two  rate  coefficients  because 
the  reactions  (Eq.  II)  depends  on  both  concentration  of  Li+  and 
oxygen  and  concentration  of  U2O2: 

JnCF=k°  (cLho2,s)  exp  -kc  (cov)  2  (co2,s)  exp 

(15) 


Dueff  =  013DLi  (7) 

Do2,eff  =  c15Do2  (8) 

Keff  =  e'5K  (9) 

rteff  =  (l-e)1-5*  (10) 


where  Du,  Dq2,  k  and  a  are  the  diffusion  coefficient  of  the  Li+  and  O2 
in  electrolyte  and  the  conductivity  of  electrolyte  and  electron  in  the 
cathode,  respectively 

2.2.  Conservation  of  charge 


Vc  =  <h  ~  fa  ~  Wfiim  -  E°  (16) 

=  jcRjilmes  (17) 

where  jc  is  local  transfer  current  density  between  electrode  and 
electrolyte  interface,  is  the  molar  concentration  of  species  i  at 
the  wall  or  surface  of  electrode,  ka  and  kc  are  the  anodic  and 
cathodic  rate  constant,  respectively,  0  is  the  symmetry  factor  equal 
to  0.5,  rjc  is  surface  or  activated  overpotential  for  reaction  at 
cathode,  A<fyym  and  Rfi/m  are  the  voltage  drop  and  the  electrical 
resistivity  across  U2O2  film  formation,  respectively,  es  is  the  volume 
fraction  of  solid  U2O2,  and  £°  is  the  theoretical  open-circuit 
potential  for  reaction. 

Whereas  the  electrochemical  reaction  at  the  lithium  metal 
anode,  it  is  given  by  the  general  Butler-Volmer  equation  as  follow 


For  the  macroscopic  model,  the  charge  conservation  between 
the  solid  phase  and  solution  phase  is  zero  and  can  be  expressed  as 
the  divergence  of  the  total  current  density  defined  by 

V-i1+V-i2=0  (11) 

During  discharge  or  charge,  the  oxidation  and  reductions 
occurring  at  the  electrode/electrolyte  interface  (charges  transfer 
reaction)  at  individual  electrode  are  represented  in  the  general 
formula  in  Eq.  (12).  The  charge  transfer  from  solid  phase  to  elec¬ 
trolyte  phase  per  unit  volume  of  electrode  (V  i2)  is  related  to  the 
average  transfer  current  density  given  in  Eq.  (13). 


ja  =  %  [e*P((1  RfnF%)  -  exp(~|r‘7?fl)]  (18) 

where  io  is  exchange  current  density  for  anode,  r}a  is  surface  or 
activated  overpotential  for  reaction  at  anode,  and  the  other 
parameter  are  as  described  above. 

The  specific  area  a  of  the  electrode/electrolyte  interface  in  Eq. 
(13)  is  decreased  by  the  morphology  and  dynamic  change  of  the 
porosity  due  to  the  L^Ch  solid  passivation  during  electrochemical 
reaction.  This  effective  local  surface  area  per  unit  volume  of  elec¬ 
trode  can  be  commonly  written  by  a  geometric  relation  [12,18] 


y^S;  Mf-me 

(12) 

V  i2  =  aj 

(13) 

This  equation  states  that  the  transfer  current  per  unit  electrode 
volume  is  equivalent  to  the  electrode  chemical  reaction  rate  where 
the  Mi  is  a  species  symbol  participating  in  the  electrochemical 
reaction,  z,-  and  s,-  are  the  charge  number  and  the  stoichiometric 
coefficient  of  the  species  i,  n  is  the  number  of  electron  transferred  in 
the  reaction,  a  is  the  specific  interfacial  area  of  the  pore  per  unit 
volume  of  the  total  electrode,  and  j  is  the  average  transfer  current 


where  eLi2o2  and  e  are  the  volume  fraction  of  solid  IJ2O2  and  initial 
electrode  porosity,  respectively.  This  empirical  equation  is  used  to 
describe  the  change  in  the  interfacial  area  for  electrochemical 
reactions  that  occur  during  discharge  because  of  the  fast  passiv¬ 
ation  of  Li202  covering  a  portion  of  the  active  sites  for  electro¬ 
chemical  reaction  over  the  carbon  surface  [19],  The  magnitude  of 
exponent  p  is  a  geometrical  factor  indicating  the  morphology 
shape  of  the  solid  peroxide  that  covers  the  active  area.  Small 
values  of  p  indicate  that  the  flat,  plate-like  precipitate  of  Li2C>2, 
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conversely,  large  values  of  p  reflects  the  needle-like  solid  which 
cover  small  active  area.  In  this  model,  the  value  of  0.5  has  been 
used. 


2.4.  Transport  through  U2O2  layer 


During  discharge  at  a  short  period  of  time,  a  thin  12202  layer 
could  be  formed  continuously,  covering  the  active  carbon  surface. 
Hence,  the  diffusion  of  Li+  and  O2  reactants  from  the  porous  media 
has  the  additional  transport  resistance  through  the  I22O2  layer.  In 
other  word,  the  active  reactants  transport  comprise  of  two  types  of 
diffusion  regimes  as  shown  in  Fig.  2.  The  first  one  is  diffusion  of  Li+ 
and  02  reactants  along  the  porous  media  and  the  other  is  their  pass 
through  the  I22O2  layer.  The  latter  transport  can  be  derived 
according  to  Fick’s  law,  which  is  proportional  to  concentration 
gradient  between  bulk  and  active  surface,  corresponding  with  the 
electrochemical  reaction  in  which  each  species  is  consumed  at  the 
active  electrode  surface  (Eq.  (14)). 


asj .  aDifilm(ci  -  cis) 
nFc  l 


(20) 


where  Dijnm  is  the  effective  diffusion  coefficient  of  species  i  across 
the  Li2C)2  film,  /  is  the  thickness  of  the  film  which  is  calculated  from 
the  volume  fraction  of  the  solid  discharge  products  given  by  Eq. 
(23).  However,  preliminary  simulated  calculation  showed  that  the 
difference  between  c,-  and  c^  is  significant  only  when  c,-  is  high.  The 
porosity  volume  change  of  the  carbon  electrode  will  be  decreased 
due  to  the  formation  of  insoluble  I22O2  covering  the  catalyst  and 
active  particles  as  described  in  Eq  (21).  Thus,  the  effective  diffu- 
sivity  for  both  Li+  and  oxygen  reactants  is  used  to  describe  how  the 
pores  are  deposited  by  the  Li2C>2  formation;  this  can  be  described  by 
the  Bruggeman  relationship. 


de 


at 


MU2 

C2FPLi2i 


(21) 


The  volume  fraction  of  the  solid  I22O2  formation  can  be  deter¬ 
mined  from  the  cathode  volume  balance  as 


eLi20  2 


(22) 


(23) 


where  Mtl-2o2  and  Pu2o2  are  the  molecular  weight  and  the  mass 
density  of  Li2C>2,  respectively,  is  the  volume  fraction  of  initial 


products 

Fig.  2.  Schematic  diagram  of  reactants  transport  inside  the  porous  electrode  and  the 
growth  of  discharge  product  on  the  surface  of  spherical  electrode  particle.  The  inset 
shows  the  enlargement  of  diffusion  transport  through  the  barrier  of  Li202  layer,  l. 


solid  phase  of  cathode  electrode  (active  carbon,  catalyst  and 
binder),  and  ro  refers  to  the  particle  radius  in  the  electrode. 

2.5.  Assumptions 

The  following  assumptions  are  used  in  this  model: 

•  The  I22O2  is  the  main  reaction  product  (Eq.  II)  and  is  only 
formed  inside  the  porous  cathode. 

•  The  electrolytes  used  in  Li— air  batteries  are  assumed  a  binary 
monovalent  electrolyte  which  consists  of  a  single  salt  in 
a  homogeneous  organic  solvent  mixture. 

•  The  Li+  diffusion  could  be  simulated  by  the  concentrated 
solution  theory.  The  pores  within  cathode  are  full  of  liquid 
phase  electrolyte. 

•  The  oxygen  is  assumed  to  dissolve  in  the  organic  electrolyte 
with  the  saturated  concentration  initially. 

•  The  convection  for  mass  transport  is  negligible  inside  the  cell. 

•  The  Li— air  cell  is  operated  in  isothermal  condition  so  that  the 
thermal  effect  is  not  considered. 

2.6.  Initial  conditions 

At  the  start  of  the  discharge  process,  initial  conditions  are 
specified  for  the  Li+  and  oxygen  concentrations,  the  porosity,  the 
specific  interfacial  area,  and  the  cell  thickness.  These  initial  values 
applied  in  the  Li— air  battery  model  are  adopted  from  literature  and 
summarised  in  Table  1.  Before  starting  the  discharge  battery,  the 
concentration  for  each  species  is  assumed  to  be  uniformly  at  all 
locations  inside  the  cell  system  and  equals  to  their  initial  concen¬ 
tration.  Because  the  thickness  of  the  porous  cathode  is  much  larger 
than  both  of  APL  and  separator,  then  this  area  is  regarded  as  the 
main  part  affecting  the  cell  performance. 

2.7.  Boundary  conditions 

From  Fig.  1,  a  schematic  view  of  the  model  cell  consists  of  four 
boundaries  and  three  regions.  The  constant  for  oxygen  concen¬ 
tration  feeding  at  the  right  side  of  the  cathode  (x  —  L )  can  be 
estimated  form  the  oxygen  solubility  (Sq,)  and  the  external 
concentration  ( c0ltext  )as  shown  in  Table  1.  At  the  current  collector 
or  the  back  side  of  the  carbon  electrode  (x  =  L),  the  current 
density  in  the  solid  phase  is  equal  to  the  applied  discharge  current 
density,  the  current  density  in  the  electrolyte  phase  equals  to  zero, 
and  the  flux  of  each  species  is  zero.  At  the  carbon  electrode/ 
separator  interface  (x  =  Lc )  the  continuous  boundary  conditions 
are  specified  for  the  fluxes  of  all  species.  The  current  density  in 
the  solid  phase  in  this  interface  becomes  zero,  and  the  current 
density  in  the  electrolyte  then  equals  to  the  applied  discharge 
current  density.  The  voltage  of  the  cell  is  calculated  by 
the  difference  between  the  electrode  potential  at  cathode  current 
collector  and  the  electrolyte  potential  at  the  anode  side, 
Vcell  =  </>i(x  —  L)  -  02(x  =  0). 

The  conservation  equations  and  the  boundary  conditions 
described  above  were  discretized  using  a  finite  element  method 
and  solved  in  one-dimensional  battery  system  by  commercial 
software  package  COMSOL  multiphysics  version  4.2a.  The  COMSOL 
software  is  designed  to  solve  a  set  of  coupled  differential  and 
algebraic  equations  and  the  battery  simulation  model  is  performed 
on  a  32  bit  Windows  platform  with  4  GB  RAM,  and  Intel  Core  2  Duo 
2.93  GHz  processor.  The  different  transport  equations  and  the 
electrochemical  reactions  were  solved  as  time  dependent.  The 
solution  were  considered  as  converged  solution  when  the  differ¬ 
ence  between  two  results  was  less  than  10  4  (relative  tolerance)  for 
all  variables. 
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Table  1 

Parameters  used  in  simulation. 


Parameter  Value 

Cell  properties 

Thickness  of  APL  5  x  10  8 

Thickness  of  separator  5  x  10  5 

Thickness  of  porous  7.5  x  10  4 

Porosity  0.73 

Specific  interfacial  3.75  x  106 

area  of  cathode 

Electrical  resistivity  50 

across  Li202  film  formation 
Electrolyte  properties 
Electrolyte  concentration  1000 

Solubility  factor  of  oxygen  0.38 

External  oxygen  9.46 

concentration 


Initial  oxygen  concentration 
at  x  =  L  (Soz*Coilext) 
Solubility  limit  of  Li202 
dissolved  in  electrolyte 
Li+  diffusion  coefficient1 
Oxygen  diffusion  coefficient 
Conductivity  of  Li+  in 
electrolyte 

Transference  number  of  Li+a 
dlnf/(dlncLi  )a 
Kinetic  parameters 
Reaction  rate  coefficient 

Reaction  rate  coefficient 
cathodic  current 
Dissolution  rate  coefficient 
of  Li202 

Exchange  current  density 
for  anode 
Symmetry  factor 
General  parameter 
Mass  density  of  Lithium 
peroxide  (Li202) 

Mass  density  of  electrolyte 
solution  (LiPF6) 

Mass  density  of  carbon 

electrode 

Operating  temperature 


KT17 


0.5 

2140 

1200 

2260 
25  x  10 

300 


a  Vary  with  concentration. 


Unit 


fim2 


kgm-3 

kgm-3 

kgm-3 

K 


Symbol  Ref. 


La  [4] 

Lc  [4] 

I  [4] 

a  [14] 

*°  [26] 

a  Calculated 

Xflm  135] 


Clip  [7] 

So,  [14] 

Co,,**  [14] 


Specific  capacity  (mAh  gc'arbon) 


Specific  capacity  (mAh  gao|id) 


Fig.  3.  Voltage-capacity  curve  on  discharge  then  charge  for  a  nonaqueous  Li-air 
battery  at  a  rate  of  0.1  mA  cm-2.  The  electrolyte  contains  1  M  LiPF6  dissolved  in  ace- 
tronitrile.  The  oxygen  solubility  factor  in  the  electrolyte  is  0.38.  The  cathode  electrode 
thickness  is  750  pm  with  porosity  of  0.73.  The  cell  cycle  is  simulated  between  2.4  and 
4.2  V  in  pure  1  atm  of  oxygen  at  operating  temperature  300  K.  (a)  The  model 
comparing  to  published  data,  (b)  The  model  comparing  to  our  group  experiment. 


3.  Results  and  discussion 

3.1.  Effect  of  applied  current  density 

The  simulated  voltage— capacity  curve  of  the  Li— air  battery 
during  discharge  and  then  charge  in  1  M  LiPF6  dissolved  in  an 
organic  solvent  between  2.4  and  4.2  V  vs.  Li/Li+  at  0.1  mA  cm-2  is 
shown  in  Fig.  3a.  It  can  be  seen  from  the  figure  that  the  simulated 
results  match  well  to  the  experimental  voltage  during  discharge  at 
1  atm  of  oxygen  [7],  During  discharge,  the  cell  potential  fell  steeply 
at  the  beginning,  from  a  voltage  of  3.4  V  to  a  plateau  at  around  2.7  V, 
and  decreased  continuously  to  2.4  V.  The  discharge  cell  potential 
around  2.5— 2.7  V  was  in  good  agreement  with  our  group  report 
previously  for  similar  battery  discharged  in  1  atm  oxygen  as  shown 
in  Fig.  3b  [20].  Although  some  of  the  parameters  applied  in  the 
model  are  different  from  experimental  data  in  Fig.  3b,  the  model 
results  and  experimental  data  showed  similarities,  suggesting  that 
the  model  is  a  promising  tool  to  identify  the  Li-air  cell  mechanisms 
and  forecast  the  cell  performance  for  design  and  scale-up.  The 
discharge  capacity  for  the  case  with  the  parameters  in  Table  1, 


based  on  the  weight  of  carbon  alone  was  about  802  mAh  g^bon.  As 
shown  in  Fig.  4,  the  discharge  process  was  terminated  by  a  rapid 
increase  of  polarisation  at  the  end  of  discharge  due  to  the  oxygen 
concentration  limitation  resulting  from  solid  Li2C>2  formation  on 
the  cathode  surface,  inhibiting  the  flow  of  reactants  (oxygen,  Li+, 
and  electrons)  to  the  active  surface.  Fig.  5a  presents  the  oxygen 
concentration  profiles  inside  the  cell  during  battery  discharging  at 
various  discharge  states.  The  oxygen  concentration  decreases  at  the 
oxygen  feed  side  of  the  cathode  and  oxygen  cannot  diffuse  further 
inside  the  electrode  because  of  the  continuous  growth  of  Li2C>2  on 
the  active  surface  of  the  porous  carbon  in  an  electrode  diminishes 
the  available  pores  for  electrolyte  resulting  in  either  increasing 
Li2C>2  insulated  film  or  pore  blocking.  Both  phenomena  lead  to 
restriction  of  oxygen  transportation  and  severely  limit  the  cell 
capacity.  To  overcome  this  problem,  the  suitable  carbon  material 
with  a  sufficiently  large  pore  size  distribution  over  the  entire 
electrochemically  active  surface  should  be  used  as  a  cathode  to 
provide  enough  space  for  the  discharge  products.  The  effect  of 
different  initial  cathode  porosities  is  presented  later.  The  variation 
in  porosity  as  a  function  of  space  and  time  for  a  discharge  current  of 
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State  of  discharge  (%) 


Fig.  4.  Change  of  positive  electrode  over-potential  during  discharge  at  a  rate 
0.1  mA  cm-2  at  different  discharge  state  (0%  =  battery  is  fully  charged). 


- 0%  Discharge  stage 

- 25%  Discharge  stage 

50%  Discharge  stage 

- 75%  Discharge  stage 

- 100%  Discharge  stage 


Dimensionless  distance 


0.1  mA  cm-2  is  shown  in  Fig.  5b.  As  the  diffusion  coefficient  of 
oxygen  is  very  low  compared  to  Li+  species  in  this  battery,  the 
porosity  falls  predominantly  at  the  oxygen  feed  side  of  the  cathode. 
As  the  extent  of  discharge  increases  the  porosity  falls  almost  to  zero 
and  thus  blocks  the  diffusion  of  oxygen  species  to  into  the  cell. 

Fig.  6  shows  the  effect  of  applied  current  density  on  the  simu¬ 
lated  discharge  curves  for  a  Li— air  battery  using  the  same  cathode 
details  as  shown  in  Fig.  3.  The  discharge  capacity  demonstrates 
a  large  decrease  from  1350  mAh  g^rbon  at  low  current  density 
0.05  mA  cm-2  to  213  mAh  g“arbon  at  a  higher  current  density 
1  mA  cm-2.  This  is  consistent  with  experiments  published  by  Read 
[7]  and  similar  trend  to  the  previous  groups  that  the  current  density 
has  significant  effects  on  the  cell  capacity  [21,22],  Moreover,  the 
discharge  voltage  plateau  is  also  lower  with  increasing  current 
density.  The  capacity  loss  at  high  discharge  rate  can  be  discussed  in 
term  of  transport  limitation  of  oxygen  diffusion  through  the 
cathode  flooded  with  electrolyte  which  cannot  maintain  the  elec¬ 
trochemical  reaction.  Thus,  oxygen  reduction  occurs  in  a  small 
region  close  to  the  cathode-current  collector  interface  as  the 
discharge  rate  increase.  Moreover,  the  rapid  porosity  reduction  due 
to  the  deposition  of  Li2C>2  on  the  surface  of  the  active  area  also 
limits  oxygen  transport  into  the  cell  and  incompletely  utilises  the 
full  capacity  of  the  electrode  porosity. 

For  battery  on  charge,  then  provided  solid  lithium  peroxide 
(eii2o2)  has  been  produced  there  will  be  a  finite  Li2C>2  concentration 
(Cn2o2 )  dissolved  in  solution  (saturation  concentration)  and  the  rate 
of  oxidation  reaction  would  be  constant  and  will  only  decrease 
when  the  solid  is  consumed  or  cL,-2o2  goes  below  the  solubility  limit. 
Therefore,  the  model  assumption  applied  on  charge  is  that  the 
dissolution  rate  (r^iss)  of  solid  peroxide  is  fast  compared  to  elec¬ 
trochemical  reaction  and  can  be  written  in  Eq.  (24).  The  dissolution 
rate  equation  is  included  for  cLi2o2  species  balance  in  Eq.  (1)  during 
battery  on  charge  and  the  driving  force  term  for  this  reaction  is 
when  the  cL j2o2  decreases  below  its  solubility  in  the  electrolyte 
(Cm*). 

rdiss  =  ^deLi202  ( Cmax  ~  CLi2o2)  (24) 


where  kd  is  the  dissolution  rate  constant  for  U2O2  during  charge 
which  assumed  very  fast,  and  cmax  is  the  solubility  limit  of  L12O2 
dissolved  in  organic  electrolyte. 


- -25%  Discharge  st 

-  -  -  50%  Discharge  st 

- 75%  Discharge  st 

- 100%  Discharge  st 


0-0-1 - 1 - , - 1 - 1 - 1 

0.0  0.2  0.4  0.6  0.8  1.0 

Dimensionless  distance 


Fig.  5.  Local  (a)  concentrations  of  oxygen  and  (b)  porosity  profile  inside  the  Li-air  cell 
during  discharge  at  a  rate  0.1  mA  cm'2  at  different  discharge  state  (0%  =  battery  is  fully 
charged).  The  electrolyte  contains  1  M  LiPF6  dissolved  in  acetronitrile.  The  oxygen 
solubility  factor  in  the  electrolyte  is  0.38.  The  cathode  electrode  thickness  is  750  pm 
with  porosity  of  0.73. 


Specific  capacity  (mAh  gj!,rbon) 

Fig.  6.  Comparison  of  the  voltage-capacity  curve  between  the  model  and  experiment 
for  a  nonaqueous  Li-air  battery  at  different  discharge  rate  (a)  0.05  mA  cm-2,  (b) 
0.1  mA  cm'2,  (c)  0.2  mA  cm'2,  (d)  0.5  mA  cm'2,  and  (e)  1.0  mA  cm'2.  The  other 
parameters  used  in  the  model  are  the  same  as  described  in  Fig.  3. 
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In  the  charge  process,  the  voltage  increased  sharply  to  reach 
a  plateau  at  about  4.0  V  as  shown  in  Fig.  3a  and  is  consistent  with 
the  published  charge  curve  at  the  same  charging  rate  with  a  slightly 
difference  [7],  The  charge  capacity  is  almost  equal  to  its  capacity 
when  discharging.  Hence,  the  model  in  which  the  published 
parameters  are  applied  can  be  used  to  describe  the  Li— air  battery 
behaviour  as  well  as  optimise  the  cell  performance  and  will  be 
compared  later  with  the  actual  experimental  data  generated  in  our 
lab.  Note  that  this  model  assumes  Li202  as  only  discharge  product 
and  does  not  consider  other  possible  parasitic  reactions,  although 
some  side  reactions,  e.g.  decomposition  of  organic  solvent,  could 
lead  to  undesired  discharge  products,  such  as  Li2C03,  LiOH,  and 
lithium  alkyl  carbonates  [23-25], 

3.2.  Effect  of  oxygen  solubility 

One  of  the  major  drawbacks  of  current  Li  oxygen/air  batteries  is 
the  low  current  densities  at  which  the  cell  operates.  The  important 
factor  that  limits  the  performance  of  Li-air  battery  is  the  solubility 
of  oxygen  in  the  electrolyte.  It  is  summarised  from  the  work  of  Read 
et  al.  that  to  improve  the  performance  of  Li— air  batteries,  one 
should  either  increase  the  diffusion  coefficient  of  oxygen  or 
increase  the  oxygen  solubility  in  the  electrolyte  to  enhance  the 
reaction  rate  of  the  cathode  so  that  the  solid  Li202  efficiently  fills  in 
the  pore  [7,10].  In  this  section,  the  solubility  of  oxygen  is  varied  in 
the  model  to  simulate  the  cell  cycling  behaviour.  In  practical,  this 
parameter  can  be  increased  by  using  solvent  with  high  capability  to 
dissolve  oxygen,  such  as  solvent  based  on  perfluorinated  solvents 
or  ether-based  electrolytes  [26],  When  the  solubility  of  oxygen  is 
increased  from  the  base  case  of  3.6  mol  m  3  to  the  higher  solubility 
of  oxygen  concentration  of  9.46  mol  m-3,  the  discharge  specific 
capacity  and  specific  energy  of  Li-air  cell  increase  as  shown  in 
Fig.  7a.  The  battery  capacity  when  discharged  at  a  low  current 
density  0.1  mA  cm-2  increased  from  802  mAh  gca1rbor|  for  the  low 
oxygen  solubility  (oxygen  concentration  3.6  mol  m~3)  to 
1400  mAh  g^al-bon  at  the  highest  oxygen  solubility  (oxygen 
concentration  9.46  mol  m-3). 

It  is  interesting  to  note  that  the  difference  in  discharge  capacity 
among  the  works  reported  previously  [7,26—28]  is  largely  a  result 
of  the  diverse  properties  of  electrolytes  to  dissolve  the  distinct 
amount  of  soluble  oxygen.  A  recent  study  by  Lu  et  al.  [26]  showed 
that  the  high  discharge  capacity  of  Li-air  battery  can  be  attributed 
to  higher  oxygen  solubility  in  the  electrolyte  with  1,2-Dimethoxy- 
ethane  (DME)  used  in  his  work  comparing  to  lower  capacity  in  the 
electrolyte  with  propylene  carbonate  (PC).  From  his  work,  the 
solubility  of  oxygen  in  the  electrolyte  with  DME  is  around 
8.76  mol  m  3  providing  the  specific  capacity  2600  mAh  at 

low  discharge  current  density  of  250  mA  Sarbon  (approximately 
0.1  mA  cm-2  in  our  simulation).  This  capacity  is  higher  than  our 
work  with  the  same  solubility  of  oxygen  concentration.  Many 
possible  explanations  could  account  for  this  but  one  factor  might  be 
that  the  diffusion  coefficient  for  oxygen  in  our  model 
(7  x  10_1°  m2  s-1)  is  lower  than  the  one  in  electrolyte  with  DME 
(4  x  10-9  m2  s-1)  and  the  transference  numbers  describing  the 
fraction  of  the  total  current  carried  by  Li+  in  a  solution  are  different 
between  the  two  electrolytes.  However,  if  these  parameters  are 
applied  in  our  model  including  the  solubility  of  oxygen  in  DME 
(8.76  mol  m-3)  at  low  discharge  current  density  of  0.1  mA  cm-2,  the 
cell  capacity  from  our  simulation  (2400  mAh  g”^,,,,)  almost  equals 
to  the  cell  using  DME  electrolyte  (the  results  is  not  shown  here). 
Therefore,  this  model  can  be  appropriately  used  to  predict  the 
behaviour  of  Li— air  battery  when  changing  detailed  parameters. 

Moreover,  it  can  be  seen  from  Fig.  7b  that  the  discharge  cell 
voltages  also  increased,  from  ca.  2.68  to  2.80  with  the  higher  dis¬ 
solved  oxygen  concentration  in  the  electrolyte.  The  cell 


energy  and  (b)  voltage-capacity  curve  on  discharge  then  charge  for  the  Li-air  battery. 
The  other  parameters  used  in  the  model  are  the  same  as  described  in  Fig.  3. 


performance  is  enhanced  because  the  oxygen  can  substantially 
diffuse  farther  inside  the  porous  structure  when  a  solvent  with 
greater  oxygen  solubility  is  used.  For  the  cell  cycling,  it  is  apparent 
from  the  same  figure  that  no  significant  difference  can  be  seen  from 
the  charge  voltage  of  Li— air  battery  at  various  oxygen  solubilities, 
i.e.  increasing  dissolved  oxygen  will  not  affect  the  cell  performance 
significantly  on  the  charging  period  as  can  be  appreciated  from  the 
kinetic  oxidation  reaction. 

3.3.  Effect  of  lithium  peroxide  solubility 

The  discharge  products  of  Li— air  battery,  Li20  and  Li202,  are  not 
very  soluble  in  a  non-aqueous  electrolyte  and  also  considered  one 
of  the  main  reasons  that  limit  the  cell  performance.  Therefore, 
introducing  some  additives  or  co-solvent  to  the  cell  electrolyte  can 
partially  enhance  the  solubility  of  the  discharge  products  and 
improve  the  cell  performance.  It  has  been  demonstrated  that 
adding  tris(pentafluorophenyl)  borane  (TPFPB)  can  substantially 
increase  the  solubility  of  Li202  in  carbonate  based  solvent  from  very 
low  amount  of  0.19  mol  m~3  to  190  mol  m~3  [29],  Consequently,  it  is 
interesting  to  use  the  model  to  simulate  the  variation  of  Li202 
solubility  on  performance.  On  addition  of  the  TPFPB  additive,  not 
only  is  the  lithium  oxides  solubility  increased  but  also  both  the  Li+ 
transference  numbers  and  the  electrolyte  conductivities  are 
enhanced  for  those  electrolytes  with  added  TPFPB  complex.  These 
effects  are  not  considered  in  our  current  model  which  is  focused 
only  the  effect  of  Li202  solubility.  From  the  data  in  Fig.  8a,  it  is 
apparent  that  increasing  the  solubility  of  Li202  in  electrolyte  does 
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not  affect  the  discharge  voltage  and  specific  capacity  of  Li-air  cell 
performance.  However,  there  is  a  significantly  improvement  of  the 
charge  cycling  cell  voltage,  i.e.  the  magnitude  of  charge  voltage 
decreased  by  approximately  0.4  V  with  higher  Li202  concentration 
dissolved  in  the  electrolyte.  The  improvement  on  charge  potential 
in  the  case  of  increasing  the  Li202  solubility  is  related  to  better 
surface  oxidation  reaction  kinetics. 

3.4.  Effect  of  kinetic  rate  constant  and  porosity 

As  the  cell  performance  depends  on  the  porous  electrode,  one 
option  is  that  the  applied  catalysts,  such  as  metal  oxide  catalysts 
and  noble  metal  catalysts,  with  carbon  active  materials  can 
reduce  the  overpotential  and  thus  increase  the  cell  efficiency.  To 
investigate  the  effect  of  different  catalysts  on  the  battery 
behaviour,  the  effect  of  cathodic  rate  constant  ( kc )  on  the 
discharge  voltage  of  the  cell  as  a  function  of  the  specific  capacity 
is  shown  in  Fig.  8b.  This  rate  coefficient  is  assumed  constant  and 
uniform  distribution  inside  the  cathode.  As  demonstrated  in  the 
results,  the  cell  voltage  plateau  gradual  increased  with  the  higher 
value  of  kc  due  to  the  reduced  overpotential  of  the  cathode 
electrode.  The  value  of  3.4  x  10-16  of  kc  demonstrated  an  onset 
voltage  of  about  2.85  V,  and  an  average  voltage  plateau  of  2.82  V 


(only  0.14  V  lower  than  estimated  equilibrium  of  2.96  V  for  Li202 
formation).  This  can  be  attributed  to  the  reasonably  high  reduc¬ 
tion  reaction  activity  at  the  carbon/catalyst  interface  with  higher 
catalyst  activity.  However,  there  was  a  relatively  small  increase  in 
specific  capacity  of  cathode  with  different  rate  constants.  This 
study  results  do  not  support  the  previous  researches  [6,27,30] 
demonstrate  that  the  various  catalysts  can  improve  both  the 
specific  capacity  and  discharge  voltage.  A  possible  explanation  for 
this  might  be  that  the  carbon  porosity  is  modified  by  the  reaction 
with  catalyst  during  the  preparation  [31],  As  a  result,  some  of  the 
closed  micro-pores  are  opened  for  the  greater  access  of  discharge 
products,  and  some  open  pores  are  widened.  However,  this  effect 
which  influences  the  cell  capacity  was  not  included  in  our  model 
with  applying  the  same  cathode  structure  (carbon  loading, 
thickness  and  porosity). 

For  the  case  of  varying  initial  cathode  porosity,  the  results  show 
that  the  battery  capacity  increases  with  the  more  space  available  to 
accommodate  discharge  products  as  presented  in  Fig.  8c.  There  is 
a  slight  increase  in  the  voltage  plateau  for  the  higher  porosity.  Note 
that  the  battery  capacity  is  normalised  to  the  weight  of  carbon 
comprised  in  the  cathode  electrode.  This  value  is  also  affected  when 
the  active  material  decreases  in  the  cell  using  high  porosity. 
Moreover,  the  porosity  is  assumed  that  there  are  no  micro-pore 
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channels  which  is  not  utilised  by  lithium  oxides  inside  the  porous 
cathode. 


3.5.  Effect  of  cathode  thickness 


Besides  concentrations  and  discharge  conditions  affecting  the 
cell  behaviour,  the  thickness  of  the  cathode  electrodes  also  plays  an 
important  role  in  the  electrochemical  performance.  Fig.  8d  shows 
the  cell  discharge  voltage  as  a  function  of  specific  capacity  for  the 
different  cathode  thickness.  Obviously,  the  thicker  electrode 
provides  the  lower  cell  performance  in  term  of  specific  capacity. 
This  can  be  attributed  to  both  the  slow  diffusion  of  oxygen  dis¬ 
solved  in  non-aqueous  electrolyte  and  the  long  residence  time 
moving  through  the  thicker  electrode.  A  number  of  experimental 
published  works  showed  that  the  scanning  electron  microscopy 
(SEM)  images  of  the  cathode  surface  on  the  airside  of  the  fully 
discharged  cell  were  almost  filled  in  the  space  by  solid  lithium 
oxides  deposition  [9,26,32],  but  this  behaviour  did  not  appear  on 
the  separator  side.  To  investigate  the  effect  of  cathode  thickness, 
the  Damkohler  number  is  introduced  as  shown  in  the  equation 
below. 


Da 


j/nF 

D02,ejrc02/L 


(25) 


The  Damkohler  number  defined  in  Eq.  (25)  expresses  a  dimen¬ 
sionless  ratio  of  current  density  of  electrochemical  reaction  rate  to 
oxygen  diffusion  rate  or  this  can  be  said  that  the  ratio  of  the 
characteristic  time  of  oxygen  diffusion  to  the  reaction  time.  A  unity 
of  Damkohler  number  means  equal  rate  of  surface  electrochemical 
reaction  and  oxygen  mass  transport  due  to  diffusion.  From  the 
Damkohler  number  for  each  cathode  thickness  in  Fig.  8d,  it  is 
indicated  that  all  the  numbers  are  greater  than  unity  and  the  higher 
value  provides  lower  specific  cell  capacity  when  discharging  Li-air 
at  the  current  density  of  0.1  mA  cm-2.  When  Damkohler  number  is 
large  (more  than  unity),  the  Li— air  performance  is  controlled  by  the 
diffusion  of  oxygen  which  has  a  low  diffusion  coefficient  in  non- 
aqueous  electrolyte  and  also  be  affected  by  the  tortuosity  of 
porous  cathode.  As  a  result,  the  cell  demonstrated  the  small  specific 
capacity  as  the  thicker  electrode  and  high  value  for  the  thinner 
electrode.  To  further  calculate  the  Damkohler  number,  Table  2 
compares  this  number  in  different  cell  thickness  and  discharge 
current  densities.  The  Damkohler  number  for  the  cathode  thickness 
of  750  pm  in  different  discharge  rate  as  shown  in  Table  2  corre¬ 
sponds  to  the  voltage— capacity  curve  in  Fig.  6  which  describes  the 
better  Li— air  performance  as  the  low  Damkohler  number  at  small 
discharge  rate.  This  study  has  shown  that,  to  meet  the  high  cell 
performance,  the  cathode  electrode  should  be  an  average  in 
thickness  around  500-750  pm  and  operates  on  the  low  current 
density.  For  the  higher  discharge  rate  or  thicker  cathode,  the  cell 
will  be  suffered  from  the  diffusion-controlled  limitation  which  can 
be  overcome  this  problem  by  using  the  electrolyte  with  high 
solubility  and  diffusivity  of  oxygen  to  decrease  the  Damkohler 
number. 


3.6.  Electrolyte  with  high  solubility  Li202  additive 

One  of  the  main  issues  during  the  operation  of  Li— air  cell  is  build 
up  of  the  solid  lithium  peroxide  depositing  and  covering  the  active 
surface  of  the  porous  cathode.  From  the  previous  section,  some 
amount  of  boron  complexes  additive,  e.g.  tris(pentafluorophenyl) 
borane  (TPFPB),  can  substantially  increase  the  solubility  of  Li202  in 
carbonate  based  solvent  within  a  limited  concentration.  Hence,  it 
would  be  interesting  to  model  the  cell  in  the  absence  of  solid  oxide 
formation  in  case  the  electrolytes  have  high  solubility  of  L12O2. 


Table  2 

Comparison  of  Damkohler  number  in  different  cell  thickness  and  discharge  current 
densities  for  oxygen  diffusion  coefficient  (D0j  =  7  x  10-lom2s-1). 


250  0.41  0.83 

375  0.62  1.24 

500  0.83  1.65 

750  1.24  2.48 

1000  1.65  3.30 

1250  2.06  4.13 

1500  2.48  4.95 


1.65 


3.30 

4.95 

8.25 

9.91 


Without  solid  formation,  the  porosity  and  active  surface  area  inside 
the  cathode  are  assumed  constant  during  discharge  operation. 
Then,  the  local  concentration  of  Li2C>2  increases  continuously  with 
time  in  the  electrolyte.  This  may  affect  the  density  and  viscosity  of 
the  electrolyte  solution  and  the  diffusivity  of  oxygen  when  high 
concentrations  U2O2  are  encountered.  The  data  for  density  of 
electrolyte  depending  on  concentration  of  Li202  are  not  known, 
therefore,  we  have  adopted  the  information  of  solution  density 
dependent  on  LiPF6  concentrations  from  [33],  The  variation  of 
electrolyte  density  is  increased  linearly  with  LiPF6  concentration 
and  the  data  is  extrapolated  to  the  desired  concentration  in  the 
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- Cell  voltage  with  solid  formation 
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Fig.  9.  (a)  Comparison  of  the  Li-air  cell  discharge  voltage  at  a  rate  0.1  mA  cm-2 
between  cell  with  and  without  solid  Li202  formation,  (b)  Local  concentrations  of 
oxygen  profile  inside  the  Li-air  cell  without  solid  Li202  formation  during  discharge  at 
a  rate  0.1  mA  cnr2  at  different  discharge  state  (0%  =  battery  is  fully  charged). 
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model.  The  volume  of  the  cell  system  is  assumed  constant  during 
simulation.  Similarly,  the  solution  viscosity  also  changes  with  salt 
concentrations  and  this  relation  is  applied  in  the  model.  Changes 
et  al.  [34],  had  investigated  the  relative  viscosity  of  concentrated 
electrolyte  solution  by  applying  the  Jones— Dole  equation  which  is 
used  to  describe  the  viscosity  of  solutions  when  the  salt  concen¬ 
tration  is  varied.  Moreover,  the  diffusivity  coefficient  of  oxygen 
decreases  with  increasing  viscosity  in  the  electrolyte.  To  describe 
this  behaviour,  the  relationship  between  viscosity  and  the  diffusion 
coefficient  of  oxygen  is  given  by  the  Stokes— Einstein  equation.  The 
results  of  cell  voltage  discharge  at  0.1  mA  cm-2  obtained  in  both  the 
model  with  and  without  of  the  solid  lithium  oxide  formations  are 
shown  in  Fig.  9a.  The  specific  capacity  predicted  by  the  model 
without  the  formation  of  solid  is  almost  equal  to  that  from  the 
model  that  includes  the  effect  of  the  porosity  change  due  to  the 
solid  formation.  However,  the  former  demonstrates  slightly  lower 
cell  voltage  compared  to  the  latter.  This  effect  may  be  from  the 
result  of  the  high  density  of  electrolyte  solution  when  highly 
soluble  salts  are  encountered. 

As  the  concentration  of  solution  increased  with  time  by 
assuming  the  high-dissolved  discharge  product  of  Li2C>2,  the  oxygen 
diffusion  rate  fell  because  of  the  high  salt  concentration  and 
increasing  solution  viscosity.  It  can  be  seen  from  the  oxygen 
transport  inside  the  porous  cathode  in  Fig.  9b  that  the  decrease  in 
oxygen  concentration  with  discharge  time  is  greater  than  the  case 
with  Li202  formation  (Fig.  5a)  and  the  diffusion  is  limited  close  to 
the  cathode/current  collector  interface  owing  to  the  high  resistance 
of  oxygen  transport  in  a  highly  viscous  solution. 

4.  Conclusions 

A  micro-macro  homogeneous  one-dimensional  model  has 
been  developed  for  the  rechargeable  Li-air  battery  using 
a  concentrated  binary  electrolyte  theory.  The  model  successfully 
predict  the  effects  of  applied  current  density,  solubility  limits  for 
both  oxygen  and  Li2C>2,  high  degree  of  Li2C>2  accumulation  and  the 
influence  of  the  cathode  structure.  This  model  considers  the  time 
dependence  and  space  dependence  of  the  battery  system  and  also 
includes  the  mass  transport  along  the  depth  of  the  cell  and  the  local 
mass  transfer  between  Li2C>2  layers  and  active  surface  morphology 
changing  with  the  Li202  growth.  The  simulated  cell  potential  for 
discharging  is  around  2.5-27  V  and  charging  at  around  4  V,  which 
agree  well  with  the  experimental  data.  The  nominal  discharge 
capacity  based  on  weight  of  carbon  alone  at  0.1  mA  cm-2  is  about 
802  mAh  grjLon,  which  is  in  agreement  with  the  experimental 
observation.  Increasing  the  solubility  limit  of  oxygen  enhances  the 
discharge  capacity  and  also  increases  the  cell  discharge  potential, 
but  does  not  affect  the  charge  potential.  Improving  the  solubility  of 
Li202  in  the  electrolyte  can  lower  the  charging  voltage  but  has  little 
effect  on  the  cell  capacity.  The  present  model  can  predict  the 
potential  and  capacity  of  the  battery  and  correlated  the  battery 
performance  to  parameters  of  reaction  species  and  cathode  struc¬ 
tures  with  a  reasonable  accuration. 
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Nomenclature 

a  specific  interfacial  area  (m2  m-3) 

c,-  concentration  of  species  i  (mol  m-3) 


c,iS  concentration  of  species  i  at  the  electrode  active  surface 
(mol  m-3) 

Cmax  solubility  limit  of  Li202  dissolved  in  organic  electrolyte 
(mol  m-3) 

Da  Damkohler  number 

Dj  diffusion  coefficient  of  species  i  (m2  s_1) 

Di  ejf  effective  diffusion  coefficient  of  species  t  (m2  s_1) 

Difiim  diffusion  coefficient  of  species  i  across  the  Li2C>2  film 
(m2  s"1) 

£  electrode  potential  of  cathode  at  any  state  (V) 

£°  electrode  potential  of  cathode  at  standard  state  (V) 

/  activity  coefficient  of  LiPF6  salt 

F  Faraday’s  constant  (96,485.34  C  mol1) 

11  current  density  in  the  electrode  phase  (A  m  2) 

12  current  density  in  the  electrolyte  phase  (A  m-2) 

I  applied  current  density  (A  m-2) 

j  interfacial  transfer  current  density  (A  m-2) 

k  reaction  rate  constant 

l  thickness  of  the  U2O2  film  formation  (m) 

La,Lc,L  thickness  of  APL,  separator,  and  porous  cathode 
respectively  (m) 
mi  mass  of  species  i  (kg) 

Mi  symbol  for  the  chemical  formula  or  molecular  weight  of 
species  i  (mol  kg-1) 

n  number  of  electrons  transferred  in  the  electrode  reaction 

Ni  molar  flux  of  species  i  (mol  m-2  s-1) 
p  surface  effect  factor 

r  particle  radius  (m) 

r,-  reaction  rate  term  in  Eq.  (1 )  that  accounts  for 

electrochemical  and  chemical  reactions  (mol  m  3  s_1) 

R  universal  gas  constant  (8.3143  J  mol-1  K-1) 

Rfiim  electrical  resistivity  across  Li202  film  formation  (Q  m2) 

Si  stoichiometric  coefficient  of  species  i  in  electrode  reaction 

So2  solubility  factor  of  O2  in  non-aqueous  electrolyte 
t  time  (s) 

t+  transference  number  of  cation  in  electrolyte 

T  temperature  (K) 

Veen  cell  voltage  (V) 

z,-  valence  of  charge  number  of  species  t 

Greek  letters 

symmetry  factor 

5  dimensionless  parameter  for  LI  /? n£/RT  (1/k  +  1/u) 

e  porosity  or  void  volume  fraction  of  porous  cathode 

es  volume  fraction  of  initial  solid  phase  of  cathode  electrode 

eii2o2  volume  fraction  of  solid  Li202  formation  in  porous 
cathode 

r)  surface  or  activated  overpotential  (V) 

k  conductivity  of  electrolyte  (S  m-1) 

Keff  effective  conductivity  of  electrolyte  (S  m  *) 
v  number  of  moles  of  ions  into  which  a  mole  of  electrolyte 

dissociates 

v+  numbers  of  moles  of  cations  produced  by  the  dissociation 

of  a  mole  of  electrolyte 

Pi  density  of  a  solid  phase  of  species  i  (kg  m-3) 
a  conductivity  of  the  electrode  (S  m-1) 

(Jeff  effective  conductivity  of  the  electrode  (S  m-1) 

0!  electric  potential  in  the  electrode  (V) 

02  electric  potential  in  the  electrolyte  (V) 

voltage  drop  across  Li202  film  formation  (V) 

V  differential  operator 

Subscripts  and  Superscripts 

1  electrode  phase 
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2  electrolyte  phase 

a  anodic 

c  cathodic 

diss  dissolution 
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